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p ' Abstract 



In this paper, we study the controhabihty properties and the Lie algebra structure of 



d ' networks of particles with spin immersed in an electro- magnetic field. We relate the Lie 

algebra structure to the properties of a graph whose nodes represent the particles and an 
edge connects two nodes if and only if the interaction between the two corresponding 
particles is active. For networks with different gyromagnetic ratios, we provide a 
^ ' necessary and sufficient condition of controllability in terms of the properties of the 

. above mentioned graph and describe the Lie algebra structure in every case. For these 

systems all the controllability notions, including the possibility of driving the evolution 
operator and/or the state, are equivalent. For general networks (with possibly equal 
gyromagnetic ratios) , we give a sufficient condition of controllability. A general form of 
interaction among the particles is assumed which includes both Ising and Heisenberg 
models as special cases. 

Assuming Heisenberg interaction we provide an analysis of low dimensional cases 
(number of particles less then or equal to three) which include necessary and sufficient 
controllability conditions as well as a study of their Lie algebra structure. This also, 
provides an example of quantum mechanical systems where controllability of the state 
is verified while controllability of the evolution operator is not. 

Keywords: Controllability of Quantum Mechanical Systems, Lie Algebra Structure, Parti- 
cles with Spin. 

AMS subject classifications. 93B05, 17B45, 17B81. 
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1 Introduction 



The controllability of multilevel quantum mechanical systems described by bilinear models 
can be investigated using results on the controllability of bilinear systems varying on Lie 
groups lTl|, [l^. In particular, general results established in can be applied to this case 



leading to the calculation of the Lie algebra generated by the Hamiltonian of the system 
and the verification of a rank condition. The determination of this Lie algebra for classes of 
quantum systems is a problem of both fundamental and practical importance in the theory 
of quantum control. In fact, it gives the set of states that can be obtained by driving the 
system opportunely and letting it evolve for an appropriate amount of time. Previous work 



in this direction, for various classes of quantum systems, was done in . 

In this paper, we analyze the Lie algebra structure and give conditions of controllabil- 
ity for a network of interacting spin ^ particles in a driving electro-magnetic field. Spin 
I particles are of great interest because they can be used as elementary pieces of informa- 
tion (quantum bits) in quantum information theory 0. These systems can be driven with 
techniques of Nuclear Magnetic Resonance ^j. A study of their controllability properties 
gives information on what state transfers can be obtained with a given physical set-up. A 



previous study on the controllability of this system was carried out in [14], |Q. Results on 
the controllability of systems of one and two spin | particles can be found in [^], |r5| . 

In the present paper we relate the Lie algebra structure of a network of spin ^ particles to 
the properties of a graph whose nodes represent the particles and whose edges represent the 
interaction between the particles. We analyze first the case of networks with particles with 
different gyromagnetic ratios. For these systems we give a necessary and sufficient condition 
of controllability in terms of connectedness of the associated graph and describe the Lie 
algebra structure in every case. It will follow from this analysis that all the controllability 
conditions are equivalent for this class of systems. In particular it is possible to drive the state 
of the system to any configuration if and only if it is possible to drive the evolution operator to 
any unitary operator. We consider then systems with possibly equal gyromagnetic ratio and 
give a sufficient condition of controllability in this case. Complete results including necessary 
and sufficient conditions of various types of controllability are obtained for low dimensional 
cases, namely for a number of particles < 3. These cases are the most common in practical 
applications. We assume here (for the case number of particles = 3) an Heisenberg model 
for the interaction between particles. In this analysis we also display an example of a model 
which is controllable in the state but not controllable in the evolution operator. 

The paper is organized as follows. In Section ^ we review general notions of controllability 
for quantum mechanical systems. We recall some results proved in 'M about the relation 



among different notions of controllability as well as some of the results of |TT|, |T2[, [113] about 



controllability of quantum systems. In Section y, we describe the general model of systems 
of n interacting spin ^ particles and define some notations used in the paper. In Section |^ 
we prove a Lemma which describes a particular subalgebra of the total Lie algebra, that we 
call the 'Control subalgebra'. This will play an important role in the following development. 
In Section |^ we study the Lie algebra structure associated to the model described in Section 
assuming that all the particles have different gyromagnetic ratios. In Section ^ we remove 
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this assumption and prove a general sufficient condition of controllability. We study low 
dimensional cases in Section |^ and give some conclusions in Section ^ 

2 Controllability of Quantum Mechanical Systems 

In many physical situations the dynamics of a multilevel quantum system can be described 
by Schrodinger equation in the form, [0], |jl8|, 

m 

\^ > = H\i^> = {A + Y^ BMt))\i^ >, (1) 

i=l 

where >[]is the state vector varying on the complex sphere 5"^^ defined as the set of n-ples 
of complex numbers Xj + iyj, j = 1, n, with Yl]=i ^| + 2/| = 1- called the Hamiltonian 
of the system. The matrices A, Bi, are in the Lie algebra of skew-Hermitian matrices 
of dimension n, u{n). If A and Bi, i = 1, m, have zero trace, they are in the Lie algebra 
of skew Hermitian matrices with zero trace SM(n)0. The functions Ui{t), i = 1,2, ...,m, are 
time varying components of electro-magnetic fields that play the role of controls. They are 
assumed to be piecewise continuous, however the considerations in the following would not 
change had we considered other classes of controls such as piecewise constant or bang bang 
controls. 

The solution of (|l|) at time t, {ipit) > with initial condition {ipQ > is given by 

|^(t)>=X(t)|^o>, (2) 
where X{t) is the solution at time t of the equation 

m 

X{t) = {A + J2BMt))X{t), (3) 

i=l 

with initial condition X(0) = /„xn- The solution X{t) varies on the Lie group of unitary 
matrices U (n) or the Lie group of special unitary matrices SU{n) if the matrices A and Bi 
in have zero trace. 

Various notions of controllability can be defined for system (|T]). In particular, we will 
consider the following three. 

• System (|Tp is said to be Operator Controllable if it is possible to drive X in to any 
value in U{n) (or SU{n)). 

^We use Dirac notation 1-0 > to denote a vector on d " of length 1, and < tp\ :— j?/" >* where * denotes 
transposed conjugate. 

■^Since trace of A and Bi, i = 1,2, m, only introduce a phase factor in the solution of (|^), and states 
that differ by a phase factor are physically indistinguishable, it is possible to transform the equation into 
an equivalent one of the same form where the matrices A and Bi, i — 1, ...,m, are skew-Hermitian and with 
zero trace, namely they are in su{n). 
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• System (P is State Controllable if it is possible to drive the state {tp > to any value on 
the complex sphere Sq~^, for any given initial condition. 

• System (0) is said to be Equivalent State Controllable if it is possible to drive the state 

> to any value on the complex sphere modulo a phase factor e**^, G R . 

From a physics point of view, equivalent state controllability is equivalent to state control- 
lability since states that differ only by a phase factor are physically indistinguishable. 

From the expression (Q) for {ip >, it is clear that state controllability is related to the 
possibility of driving X to a subset of SU{n) or U{n) which is transitive on the complex 
sphere. Transitivity of transformation groups on spheres was studied in 0, |T6[, [[l^, |]20[ 
and the necessary connections for application to quantum mechanical systems where made 
in |T]]. In the following theorem, we summarize some of the results of [|T| that will be used in 



the following. Part 2) of the Theorem was proved in [jTT|, [|T^, |T^. Here and in the following 



we will denote by C the Lie algebra generated hj A, Bi, . . . , Bm in (|l|). 
Theorem 1 

1. A quantum mechanical system ^ is state controllable if and only if it is equivalent 
state controllable. Both these conditions are implied by operator controllability. 

2. The system is operator controllable if and only if the Lie algebra C generated by the 
matrices A, Bi, Bm is u{n) or su{n). 

3. The system is state controllable if and only if C is su{n) or u{n), or, in the case of n 
even, isomorphic to sp{^)^. 

4. Consider the n x n matrix with i in the position (1, 1) and zero everywhere else. Call 
this matrix D. Let V be the subalgebra of C of matrices that commute with D. Then, 
the system is state controllable if and only if dim. C — dimP = 2n — 2. 

5. Assume n even. There is no subalgebra of su{n) which contains properly any subalgebra 
isomorphic to sp{'^) other than su{n) itself. 

Because of the equivalence between state controllability and equivalent state controllability, 
in the sequel we will only refer to the two notions of state controllability and operator 
controllability. Controllability notions in a density matrix description of quantum dynamics 
were considered in [111. 



■^Recall the Lie algebra of symplectic matrices sp{k) is the Lie algebra of matrices X in su{2k) satisfying 

h 



XJ + JX^ = 0, with J given by J = 
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3 Model of interacting spin | particles 

Prom this point on, we will denote by n (which in the previous section denoted the dimension 
of a general quantum system) the number of spin | particles in a network. The state 
dimension of this system is 2". 

To define the model we will study, we first need to recall some definitions. The following 



three matrices in su{2) are called Pauli matrices (see e.g. [|T9 



1 /O 1\ 1 /O -i\ 1 /I \ 

^^^=2ll Oj' ^^^=2^ oj' 2^ -J- 

The Pauli matrices satisfy the fundamental commutation relations 

[o-^., ay] = ia^; [dy, a^] = ia^; [a^, a^] = iay. (5) 

It is known that the matrices iax, iuy, iaz form a basis in su{2). Moreover, the set of matrices 
i{cri ® (8) ■ ■ ■ C?) cr„), where aj, j = 1, ...n, is equal to one of the Pauh matrices or the 2x2 
identity /2x2, without z(/2x2 ® -^2x2 (S> ■ ■ ■ ® /2X2), form a basis in sm(2"). (Here ® indicates 
the Kronecker product for matrices.) 

In the following, we will use the notation 1^^. for the Kronecker product 



' kx 



cTi ® 0-2 ® ■ • ■ ® a„, (6) 



where all the the elements aj, j = 1, ...,n are equal to the 2x2 identity matrix, except the 
k—th element which is equal to ax- More in general, we will use the notation Ikih,k2h,---,krir^ 
with 1 < ki < k2 < ■■■ < kr < n and Ij = x,y or z, j = 1, r, for a Kronecker product of the 
form (^ where all the aj are equal to the identity /2X2 except the ones in the kj—th positions 
which are equal to the Pauli matrices ai.. The matrices so defined (excluding the identity 
matrix) span su(2"). Some elementary properties of the commutators of the matrices just 
defined that will be used in the following are collected in Appendix A. 

The Hamiltonian of the system of n interacting spin ^ particles in a driving electro- 
magnetic field is given in the form [^: 

H = Ho + Hi. (7) 
Here Hq, which denotes the internal (or unperturbed) Hamiltonian, is given by 

n 

Ho '■= y^X^klIkx,lx + Nkllky,ly + Pkllkz,lz), (8) 
k<l 

where M^i, N^i, P^i are the coupling constants between particle k and particle /. This general 
model of the interaction between different particles includes as special cases both the Ising 
and the Heisenberg model ([0, pg. 46). The term Hj, Control Hamiltonian, is given by 

n n n 

Hi := CZl lklkx)ux{t) + (XI lkhy)uy{t) + (^ Jkhz)uz(t), (9) 

k=l k=l k=l 
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where m^, Uy and Uz are the x, y and z components of the electro-magnetic field and 7^-, 
j = 1, 2, is the gyromagnetic ratio of the j-th particle. In general, we assume that we are 
able to vary all the three components of the magnetic field for control (cfr. Remark |5.2| ). 
Schrodinger equation for the evolution matrix X has the form, 

X = AX + B^Xu^ + ByXuy + BzXuz, (10) 

with 

n 

A := —i ^ {Mkllkx,lx + Nkllky,ly + Pklhz,lz), 
k<l,k,l=l 

and 

n 

Bv ■= -iC^lklkv), with V = x, y, or z. 

k=l 

It is clear that the controllability properties of this class of systems only depends on the 
parameters Mki, Nki, Pki and 7^. Our goal in the next sections is to characterize the structure 
of the Lie algebra generated by A and B^, By, Bz, C, in terms of these parameters. The 
network of spin particles can be represented by a graph whose nodes represent the particles 
and are labeled by their gyromagnetic ratios and an edge connects the nodes corresponding 
to particles k and / if and only if at least one of the coupling constants Mm, Nm, Pm is 
different from zero. In this case, the edge is labeled by the triple {Mki, Nki, Pki}- It is our 
goal, in the next sections, to relate the properties of the Lie algebra £, generated by A, B^, 
By and Bz to the properties of this graph. In the following, we denote this graph by QV. 

We define an ordering on the n particles so that the first ni have the same gyromagnetic 
ratio 7i, the next n2 particles all have gyromagnetic ratio 72, with 72 7^ 71, and so on up 
to the r— th set of Ur particles with gyromagnetic ratio 7^, with 7^- 7^ 7^^. when j ^ k and 
rii + n2 + ris + ■ ■ ■ + rir = n. We shall denote the first set of particles by S^, the second one 
by 5*2, and so on up to the r— th, 5°. We also define, for j = 1,2, ...,r, v = x, y, z 

Ijv '■= X! ^^1-^^ (ll) 

and we have 

r 

i=i 

For a given system, we shall call the Control Subalgebra of C, the subalgebra generated by 
the matrices B^, By and Bz- We shall denote the control subalgebra by B. 

4 Characterization of the Control Subalgebra 

The following lemma shows that the control subalgebra S of a spin system is the direct sum 
of r subalgebras isomorphic to su{2). 
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Lemma 4.1 Assume we are given a model as in (|T0D, and let 71, ... ,7^ be the different 
values for the gyromagnetic ratios. Assume that to each value jj correspond rij particles in 
the set 5'°, j = 1, . . . , r, then the matrices B^, By and i?^ generate the following Lie algebra: 



B = B,®By(BB,, (12) 



with: 



B^ = spanj=i^,„^r{iljx}, (13) 

By = spanj=i^_^r{iljy}, (14) 

B^ = spanj=i^,„^r{iljz}- (15) 

Moreover, we have: 

[B,, By] = B,, [By, B,] = B^, [B,, = By. (16) 

Proof. First, notice that Ij(x,y,z) satisfy the commutation relations 

where we used the Kronecker symbol 5jk. We proceed by induction on r > 1. If r = 1, then 
we have, for v G {s, y z}: 

By 

thus (p^-(p^) follow immediately from the basic commutation relations (0). 
To prove the inductive step, we first show, again by induction on r > 1 that: 

[B.,By] = -2E,^=i7|4> 

[By,B,] = -2E-=i7|/,x, (18) 
[B,,B,] = -2E-=i7|4- 

We will prove only the first of the previous equalities, since the other ones may be obtained 
in the same way. If r = 1, then 

[B^,By] = --fl[ii^jiy] = -i-ffiiz, 

where to get the last equality we have used (|l3). Now let r > 1: 

r r 

[B^,By] = -[J2^j^jx,J2^j^Jy] = 

(r— 1 -i — 1 r— 1 r— 1 

Ij^jxi ^ lj^jy\ + i^' 7r-^ry] + ^[7r-^ra) Ij^jy] + [pirlrxi Ir^ry] 

j=l j=l j=l j=l 

By the inductive assumption, we have: 

i — 1 r—l i — 1 

E E = ^ E ^j^jz- (19) 

j=i j=i j=i 
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Using (p!?!), we obtain, for j < r, 



and 



Now putting together equations (p^D, ( pOD and (|21|), we get: 

as desired. Thus, we have proved (|TBp. 

Now notice that, for example, [By, B^] has the same form as B^ except that the 7j's have 
been replaced by 7^, therefore, using the same arguments as above one may show that: 

[[By,B,],By] = ~zj2^p,^. (22) 

More in general, considering the Lie bracket between := —iY7j=il^Ijxi and Gy := 
—iJ2^j=i'jjIjy, we get 5* := —iYJj=il^^^Ijz- Proceeding this way, we obtain all the ma- 
trices 

^jll]hx. (23) 
i=i 

^jZihy^ (24) 

and 

^E7]4, (25) 

/ = 1, ...,r. The matrices in ( p3D form a basis in since the Ij^ do and the linear trans- 
formation in (^) is nonsingular. In fact, the corresponding determinant is a Vandermonde 
determinant which is different from zero because all the 7j's are different from each other. 



The same is true for the elements in (^) and (25) which form a basis in By and Bz-, respec 



tively. Finally, the commutation relations (|T^ follow immediately from (|r^). □ 



Notice that it follows from ([T^) and (^ that the subalgebras spanned by Ij[x,y,z) are 
each isomorphic to su{2) and they commute with each other. For a given j, the Lie group 
corresponding to span{Ij(^x,y,z)} is given by rij copies of SU{2) (where rij denotes denotes 
the number of particles with gyromagnetic ratio 7^) 0. Therefore it is isomorphic to S0{3) 
or SU{2) according to whether nj is even or odd, respectively. 

"^This is the Lie group of matrices of the form /i eg) L ® /2, where the identity matrix Ii has dimension 
2ni+ - nj-i^ the identity matrix I2 has dimension 2"~"i~"-2 "3 a,nd L has dimension 2"^ and is has the form 
F F ® ■ ■ ■ ^ F, with F G SU{2) and the Kronecker product having nj factors. 
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5 Lie Algebra Structure and Controllability with Dif- 
ferent Gyromagnetic Ratios 

In this section, we shall assume that the gyromagnetic ratios 71,. ..,7^ are all different. 



Therefore we have r = n and, from Lemma [4.1|, we have that the control subalgebra B is 



the span of the ilj(x,y,z), j = 1, ■■■,n. We shall give a necessary and sufficient condition of 
controllability and describe the nature of the Lie algebra C, in terms of the properties of the 
graph QV. This graph will, in general, have a number s of connected components. We first 
describe the situation when s = 1 and then generalize to the case of arbitrary s. 



Theorem 2 Assume we are given a model as in j^TSj), where the values 'jj, j = 1, . . . ,n of 
the gyromagnetic ratios are all different. If the graph is connected, then 

C = su(2"). (26) 

As a consequence the system is operator and state controllable (see Theorem^. 

Proof. We show that all the matrices of the form ilkihMh^-^kmim, can be obtained as repeated 
commutators of A, B^, By., Bz, for every 1 < m < n. Lemma [4.1| gives the result for m = 1. 
We first prove that this is true for m = 2 as well, and then proceed by induction on m. 
If m = 2, we want to show that we can obtain all the matrices of the form ilkv,iw, k < I, 
V, w E {x, y, z}. From our assumption on the connectedness of QV, there exists a path 
joining the node representing the k — th particle and the node representing the /— th particle. 
Let us denote by p the length of this path, namely the number of edges between k and /. We 
proceed by induction on p. If p = 1, then at least one among M^i, N^i and Pm is different 
from zero. If Pki 7^ 0, we have: 

[A, ill,] = i Y.(-Nhlhy,lz + Phlhz^ly) + Y.(-Nlhllz,hy + Plhlly,hz)^ (27) 
\h<l h>l J 

and 

[[A,ili,],-ilky] = -iPkihxiy (28) 

Since Pki 7^ 0, from the matrix —iPkilkxiy, using (repeated) Lie brackets with elements il^f 
and/or ilif, with /, /' G {x, ?/, z} one can obtain all of the elements of the form ilkv,iw, with 
V, w E {x, y, z}. If Pki = 0, but A''^^ 7^ 0, the same can be proved by taking the commutator 
with ilix first and then the commutator with ilkz and analogously, if Nki = Pki = 0, by taking 
the commutator with iliy first and then with ilkz- Now, assume it is possible to obtain every 
ilkv,iw for every k < I whose distance is < p — 1. Let k and / have a path with distance p 
and let / represent a particle/node in between k and I in the path. Let us also assume just 
for notational convenience that k < I < I. From the inductive assumption, we know that 
ilkvjw and Hff^fi can be obtained for every v, w, /, /' G {x, y, z}. We need to show that 
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we can also obtain every ilkg,iq for every g, q E {x, y, z}. Using equation (|66D in Appendix 

A, we get 

and 



ihx,lz,ly^ihz,lx] — -ilkx,lz, (30) 



1 

4 

where we have used the folfowing property of the Pauh matrices 

= = = 4^2x2- (31) 

As before, we can now take repeated Lie brackets of the matrix obtained in (^) with matrices 
of the form il^f and/or ilif, with /, /' G {x, y, z}, to obtain all of the matrices ilkv,iw, for 
V, w E {x, y, z}. This concludes the proof that every Kronecker product with two matrices 
different from the identity can be obtained, namely m = 2 in the above notations. 

We now show that every matrix iIkiviMv2,---,k„^vm can be obtained. Consider the Lie 
bracket 

\~'^IklV\,k2V2,...,km-lX-!T'Ikm-\y,k,nVm\ ~ '^■^k\V\,k2V2,...,km-lZ,kmVm- (32) 

Both elements —iIkrui,k2U2,.--,km-ix and ilkm-iy,kmvm available because of the inductive 
assumption. If Vm-i = we have concluded otherwise, the Lie bracket with the matrix 
ilkm-ix or ilk^_iy leads to the desired result. This conclude the proof of the Theorem. □ 

In the general situation, assume that QV has s connected components and denote by rij 
the number of nodes in the j— th component. Set up an ordering of the particles so that 
the first rii are in the first connected component of the graph, the ones from rii + 1 up to 
rii + n2 are in the second component and so on. We have ni+n2 + ■ ■ -ng = n. The following 
Theorem describes the structure of the Lie algebra C in the general case assuming to have 
different gyromagnetic ratios 7^, i = 1, 2, n. 



Theorem 3 Assume we are given a model as in ^Tuj), where the values 7^, j = 1, . . . 
of the gyromagnetic ratios are all different. Moreover, assume that the graph QV has s 
connected components (as described above), then 

£ = 5i © 52 ® ■ ■ ■ © 5,, (33) 

where each Sj, j = 1,2, ...,s is the subalgebra spanned by the matrices 

HklVl,k2V2,...,krVr : (34) 

with 

ni + n2 + ■ ■ ■ + Uj^i < ki < k2 < ■ ■ ■ < K < ni + n2 + ■ ■ ■ + Uj . (35) 
Proof. First notice that, from equation ([6^ ) in Appendix A, it follows immediately: 

[Sj,Sk] = 0, iijj^k. (36) 
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Since the values 7^ are all different, from Lemma ^TT] we have that all the elements of the 
form ilkv, k = 1, . . . ,n, V & {x, y, z}, are in C. We can write the matrix A as 



l<k<l<nx{MklIkxlx + Nkllkyly + Pkllkzlz) + 

{Mklhxlx + Nkllkyly + Pklhzlz)+ (37) 
■ ■ ■ + J2ni+n2 + ---ns^i<k<l<n{MklIkxlx + Nkllkyly + Pkl^kzlz), 

using the fact that M^/ = A^^/ = P^; = if and Z are in two different connected components. 
Taking the Lie brackets with elements ilkv, v G {x, y, z}, with rii + n2 + ■ ■ -"^j-i < k < 
ni + n2 + ■ ■ -fij (here if j = 1, we put no = 0), one may show, as in the proof of Theorem 
0, that it is possible to obtain all the elements in Sj, j = 1,2, ...,s. Moreover from (p^), it 
follows that these and their linear combinations are the only matrices that can be generated 
by A, B^, By, B,. □ 

Notice that, in the above situation, one may think of the spin system as a parallel 
connection of s spin systems of dimension rij, j = 1 . . . , s, controlled in parallel by the same 
control. The solution of (0) has the form 

X{t) = <l>i(t)<l>2(t) • • • $.(t), (38) 
where $j(t) is the solution of (|iy) with 

A = -i {Mhkhx,kx ~l~ ■Nhk-^hy,ky + Phklhz,kz)i (39) 

nj_i<fe<fe<nj 

and 

Bv = -i Y Ikhv, V e {x, y, z}. (40) 

k=nj-i+l 

The controls are the same for every subsystem and the matrices $j in ( |38D commute due 
to (pGl). The set of states that can be obtained with an appropriate control for system 
( PUf ) is given by the Lie group corresponding to the Lie algebra £ namely, in this case, 
5t/(2"i) (g) 5f/(2"2) ® ■ • • O 5f/(2"0. 



Remark 5.1 It is important to notice, and it will be used later in the next Section, that, 
in Theorems ^ and |^, the assumption of different gyromagnetic ratios is used only to derive 
that the Lie algebra spanned by ilj(x,y,z) is a subalgebra of C Thus both statements of 
Theorems ^ and ^ remain true if, instead of assuming 7^ 7^ 7^ for all i ^ j, we assume 
spanj=i^,,,^n{'ilj{x,y,z)} ^ This fact will be used in the following Section. 



In the following Theorem, we answer the question of state controllability for spin systems 
with different gyromagnetic ratio. It follows from Theorem |l] that, if £ = sm(2"'), the set 
of states reachable for system (|I^) is SU{2"-) and therefore the system is both operator 
controllable and state controllable in this case. If £ 7^ S'u(2"'),we have seen that the set of 
states reachable for (|ToD is SU{2'^^) ® 5f/(2"2) ® ■ ■ ■ ® 5[/(2"=). To see that the system is 
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not state controllable, notice that the corresponding Lie algebra C is not simple (since each 
of the subalgebras isomorphic to su(2"^ ) is actually an ideal in C) and therefore it cannot be 
isomorphic to sp(2"~^) as in Theorem 1, part (3). A more direct and geometric proof of the 
fact that S'?7(2"'i) ® S'?7(2"'2) ® ■ ■ ■SU{2^'') is not transitive on the complex sphere is to reason 
as follows. Assume for simplicity s = 2 and Vi and V2 two subspaces, of dimension 2"^ and 
2"2 such that the underlying subspace of the overall system is V^i ® V2. Every 'unentangled' 
state, namely a state of the form |t>i > ®|f2 >, with vectors |f 1 >G Vi and \v2 >G V2 can only 
be transformed into another unentangled vector {A (g> B){\vi > (X>|f2 >) = ^l^i > ^B\v2 > 
and there is no possibility of transforming \vi > (S>|f2 > into an entangled vector namely a 
vector that cannot be written as the tensor product of two vectors from Vi and V2. On the 
other hand, entangled states always exist for a pair of non trivial vector spaces Vi and V2 
(for example, if \ej >, j = 1, mi, is a basis of Vi and \fk>,k = l, 1712 is a basis of V2, 
so that \ej > \fk > is a basis of Vi ® V2, consider > |/i > +-^\^mi > \fm2 >•) We 

summarize the results in this section with the following theorem. 

Theorem 4 Consider a system of n-spins with different gyromagnetic ratios given by the 
model ^Tdj). For this system all the controllability notions are equivalent and they are verified 
if and only if the associated graph QV is connected. 



Remark 5.2 In many physical implementations of the control of spin | particles, the z 



component of the control is held constant, 
occur in the proof of Lemma |0. 



The only changes in the previous treatment 
In fact, for this case, one does not have the matrix B^. 
However, by using the first one of equations (0), one obtains —iY7j=il'^Ijz G -B. Then, 



using this matrix in place of B^, one gets all the matrices in (p3|). 



Vs in 



4]), and even /'s in (^51). If we assume |7i| 7^ I72I 7^ • • • |7r 



(p5D , with only odd 
, the result remains 



unchanged. In fact, the determinant of the matrix referred to at the end of the proof of 
Lemma 4J., is still a non zero Vandermonde determinant. The drift matrix A is modified by 
adding a term —iJ2]=i'ljIjzUz, with Uz constant but this does not modify the resulting Lie 



algebra C, since —i Yl]=i IjljzUz belongs to the control subalgebra. 



6 Systems with Possibly Equal Gyromagnetic Ratios 

In this section we analyze the graph QV for networks of spins with possibly equal gyro- 
magnetic ratios and give a sufficient condition of operator controllability for these systems 
in terms of the properties of this graph. It will follow from the examples in the next sec- 
tion that the equivalence between state controllability and operator controllability, proved in 
Theorem ^ for systems with different gyromagnetic ratios, does not always hold if we allow 
two particles to have the same gyromagnetic ratio. 

In the following we describe an algorithm on the graph to conclude operator con- 
trollability. The main idea and the physical interpretation go as follows. When all the gyro- 
magnetic ratios of the particles are different they 'react' in a different way to the common 
electro-magnetic field and this 'asymmetry' along with connectedness of the spin network 
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allows us to control all the particles at the same time. However, even if two particles have 
equal gyromagnetic ratios they might interact in different ways with a third particle which 
has gyromagnetic ratio different from the two, and this will break once again the symmetry 
and give controllability. 

Let us divide the particles into r sets 5*°, 5° as it was done in Section 3 and assume that 
at least one set is a singleton, namely, there exists at least one particle which has different 
7 from all the others. Consider a set S containing all the singleton nodes. Assuming that 
there are m of them, let the sets S'5',...,S'°_^ be of cardinality > 2. Now we illustrate a 
'disintegration' procedure to divide these sets further. 

Algorithm 1 

1. Let C be a collection of sets. Set C := S^, S2, S^_^. 

2. For each set S in C, consider a particle / in S such that for at least two particles in k 
and j in S 

{\M,i\, \N,j\, \P,j\} ^ {\M^j\, \N^j\, |P^,-|}. (41) 

If there is no element in S and no set in C having this property STOP. Divide the set 
S into subsets of particles that have the same value for {|M;j|-|, lA*"^;-], |-Pm|}- 

3. Consider the sets obtained in Step 2. Put the elements that are in singleton sets in S. 
If all the elements are in S, STOP. 

4. Replace the collection C with the remaining non singleton sets and go back to Step 2. 

We have the following theorem. 

Theorem 5 // Algorithm 1 ends with all the particles in the set S and is connected, 
then the Lie algebra C associated to the spin | particles system, with n particles, is su{2^). 
As a consequence the system is operator controllable. More in general, if Algorithm 1 ends 
with all the particles in the set S and QV has s connected components of cardinality ui, n2, 
ris, C is given by (See Theorem^. 

Proof. From Remark |5.1| , all we have to show is that, in the given situation, the Lie algebra 
spanj=i^^^^^n{'ilj{x,y,z)} is a subalgebra of C. Rewrite the drift matrix A as 

A = —i ^ {Mkllkx,lx + Nkllky,ly + Pklhz,lz) 

—i X! {Mkllkx,lx + Nkllky,ly + Pklhz,lz) 

— i X! {Mkllkx,lx + Nkllky,ly + Pklhz,lz) 

—i X! {Mkllkx,lx + Nkllky,ly + Pklhz,lz)- (42) 

k<l,k^S° ,l€S° 

' r — m' r-~m 
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From Lemma |4.1| , the matrices ilj^, v G {x, y z} and j = 1, 2, r, where r is the number of 
sets Sj, are available to generate the Lie algebra C. In particular, since we have assumed that 

the last m sets are singletons, the matrices i/^, v G {x, y, z}, I = ni + n2-\ h nr-m+, n 

are G C Now, assume that in the set S'°_^ there are two elements j and k such that 
condition ( ^T]) is verified for some / G S and assume, for the sake of concreteness, that the 
inequality is verified for the P coefficient (minor changes are needed in the other cases). By 
taking the Lie bracket of A with il(r-m)x, the first term gives zero, since it does not involve 
any term in the set (see (|6^, in Appendix A and the definition of the J's in (p!TD). The 

Lie bracket of the second term with il(r-m)x gives a matrix which is a linear combination 
of matrices of the form ilkv,pw, k,p E S^_^ and v,w E {x,y, z}. We call this matrix Kr-m- 
Thus, we have 

{A, il(^r—m)x^ -^r— m~l~ 

^ {T.k<l,keSO_^,l(^S'^_^i-NklIkz,ly + Pklhy,lz) + Y^k<l,kfS°^^,leSO^^i-NklIky,lz + Pklhz,ly)) ■ 

(43) 

By taking the Lie bracket of ( ^3|) with i/^^, and using Properties 1 and 2 in the Appendix 
A, we obtain 

[[A,ii(^r-m)x],iPly]=i Yl Pkihyjx- (44) 

From this matrix, by taking Lie brackets with il(^r-m)v and/or ilj^, v G {x, y z}, it is possible 
to obtain all the matrices of the form ( PI) with all the possible combinations of x, y and z 
in place of y and x respectively. 

Using ( p^ ) in Appendix A, it is not difficult to see that 

X! Pkl^IkyJz-,'^' X! Pkl}IkxJz] = X! ^kl^kz- (45) 

k&S2_^ k^S2_^ k^S2_^ 



r — m 



By taking the Lie bracket of this with —i I]fce5" Pkf^kxjz^ obtain i I]fces" Pkfhyjz 
and repeating the calculation as in (^5|) , we obtain 

Pkfhyjz,i Y Pkiihxjz] = \i Y Pkjhz- (46) 

Continuing this way, it is possible to obtain all the matrices of the form 

' E Plllkz. p = 0,l,2,..., (47) 

and, with minor changes in the choice of the Lie brackets, we can obtain 



I E Pllhx, ^ E Plpky, p = 0,1,2,.... (4^ 
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Now consider for example, the matrices iJ^kGS^ Pjikz and assume, without loss of gen- 

r — m 1^'' 

erality that the elements P^flkz are arranged so that elements that have the same value for 
Pki appear one after the other in the sum. The associated determinant is (cfr. the proof of 
Lemma |4.1| ) a Vandermonde determinant and therefore by appropriate linear combinations 
we can obtain all the matrices of the form X^feeT hz where T is a generic subset of S^_^ such 
that all the values of \Pki\ are the same, for all the k & T. In particular, if T contains a 
single element then we place that element in the set of singletons S. The other subsets of 
Sr-m arranged in new sets. It is clear that we can repeat this procedure for the other 
sets Si, S2, Sr_jj^_i, and then for the subsets obtained, as described in Algorithm 1. If 
the procedure ends with all the elements in S then we have that spanj=i^,,,^n{'ilj{x,y,z)} is in 
C and the Theorem follows from Remark □ 



7 Low dimensional systems 

Results on the controllability of spin systems in the cases of n = 1 and n = 2 particles can 
be found in 0, |^ and In this section we consider the model ([T0|) assuming Heisenberg 
type of interaction namely 

Mki = Nu = Pki-= Jkh (49) 

for every pair of particles k and /. For this model, in the case n = 2, the only noncontrollable 
case, C 7^ su{2), is when n = 2 and the two particles have the same gyromagnetic ratio. In 
this situation, we have 

C = span{A} ® span {i {a^ ® I2x2 + hx2 ® <^v) , ve{x,y,z}}, (50) 

and the matrix A commutes with all the matrices in C. The Lie algebra C is isomorphic to 
ui2). 

We treat now completely the case of n = 3 interacting spin | particles. If the three 
particles have all different gyromagnetic ratios, then we are in the situation treated in Section 
^. There are two more possibilities: 

(a) all the three gyromagnetic ratios all equal (i.e. 71 = 72 = 73), 

(b) two gyromagnetic ratios are equal and the third is different (i.e. 71 = 72 and 71 7^ 73, 
according to the notations in Section H, we have 5° = {1, 2} and 5'° = {3}). 

• case (a) 

This case is particularly simple. In fact, we have: 

C = span{A} Q) span{iliv, v&{x,y,z}}, (51) 

with 



span{A}, span{ilix, ihy, ihz} 



0. 
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The Lie algebra C is isomorphic to m(2) and the model is neither operator controllable nor 
state controllable from Theorem 

• case (b) 

This situation is more involved and it gives rise to interesting examples. First recall that, 
from Lemma [4.1| , we get, for v = x,y, z, 

= Span{ -Z (cr^ O /2x2 + ^2x2 ® CTt;) ® ^2x2, -«(/2x2 ® ^2x2 ® Cr«)}, /c.<.x 

B = B,®By®B,. ^^^> 
To deal with this case, we need to consider three subcases: 

(i) iJisI ^ IJ23I, 

(ii) Ji3 = J23, 

(iii) Ji3 = -Jaa- 

For the case (i) we can apply Theorem || and conclude that, if the associated graph is 
connected then C = su{8) and the system is operator controllable. For the case (ii), the 
model will turn out to be neither operator controllable nor state controllable. Finally, in the 
case (iii), the controllability properties of the model will depend on the coefficient Ju- In 
fact the system will be operator controllable (i.e. C = su{8)) if J12 7^ 0, while, if J12 = 0, 
then the system will be state controllable but not operator controllable (so, from Theorem 
0, in this case £ is isomorphic to sp(4)). 

• case (ii): J13 = J23 

From a physical point of view, in this case the particles one and two feel the same magnetic 
field and have the same interaction with the third particle, therefore it is not possible to 
manipulate separately these two particles. This internal symmetry of the system results in 
lack of controllability both for the evolution operator and the state. If J13 = J23 = 0, we 
have: 

• if J12 = 0, then C = B, 

• if J12 7^ 0, then C = span{A} © B and the matrix A commutes with all the matrices 
in C. 

Now we consider the case J13 = J23 7^ 0. We first define an operation of 'symmetrization' p 
on the matrices in m(4), as follows: 

ip (cTi (g) 0-2) = (cTi (g) 0-2 + cr2 ® cTi) , (53) 

with o"i,(T2 G {l2x2,o'x,o'y,o'z}, and we extend p to all of the matrices of m(4) by linearity. 
Let: 

J'p = {Xe m(4) I p{X) = X} . (54) 
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Notice that: 

Xi., X2 E J^p ^ X1X2 ^ Tp. 
For sake of completeness, we include a proof in Appendix B. Let: 



n = < 



H = F' 



pi 



C^i G {/2x2, (^z) 



> . 



First, we have: 

ccn. 

To see this, recall that C is generated by: 

A = -iJi2 ((Tx- ® CTnc ® hx2 + cTy^ay® /2x2 + cT^ ® 0"^ ® /2X2) - i-hz 

((cTx ® hx2 + ^2x2 ® ® 0-x' + (0-3; ® ^2x2 + ^2x2 ® O"?;) ® (T^; + (o"^ (g) /2x2 + ^2x2 ® CT^) 



(55) 



(56) 



(57) 



and by the matrices in B (see equation (^21) ). Thus £ C 7i follows from the fact that both 
A and B are in 7i, and that 7i is a Lie algebra because of (|55|). Now we have: 

(i) if Ji2 7^ 0, then C = Ti, and it has dimension 39; 

(ii) if Ji2 = 0, then CcTi, where the inclusion is strict and it has dimension 38. 

The proof of both the previous statements (i) and (ii) follows from the analysis of the Lie 
algebra structure for this model, given in the Appendix C. In both cases C is not su{8), thus 
the model is not operator controllable. Moreover, by looking at the two possible dimensions 
of £, the model can not be state controllable either. In fact to have state controllability we 
would need, see Theorem 0, £ = su{8) or C isomorphic to sp{4), which has dimension 36. 

• case (iii): J13 = -J23 7^ 

This case is interesting because it provides a physical example of a system which is state 
controllable but not operator controllable. It also shows that for spin systems with some 
gyromagnetic ratios possibly equal to each other the two notions of controllability do not 
coincide (cfr. Theorem^. 



Consider the following vector spaces of matrices 

M := span{ili^^3^ - il2v,3w, v,w e {x, y, z}}, 
C := span{iliy + il2v, ihw, v,w e {x, y,z}}, 
M := span{iliy^2w,3p + ihw,2v,3p, v,w,p E {x,y,z}}, 
TZ := span{ili^^2w - ihw,2v, v j^w,v,w e {x, y, z}}. 



(58) 
(59) 
(60) 
(61) 

It can be seen by verifying the commutation relations among these vector spaces that A := 
A^©C©A/'©7^isa subalgebra. Moreover, using the test in part 4 of Theorem |I], it can 
be shown that this Lie algebra is isomorphic to sp(4). It is interesting to notice that the 
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decomposition A := A4 (B C (B (B TZ is underlying a Cartan decomposition of sp(4) |jTO| 
since the following inclusions among the above vector spaces hold: 

[c®^f,c®^^]cc®^^, [c®^^,M®n]cM®n, [M®n,M®n]cc®^^. (62) 



To see that ^ is a subalgebra of C notice that Lemma [4.1| gives a basis for C. By taking the 
Lie bracket of A with /3a; G B and then of the resulting matrix with I^z G B, we obtain a 
matrix proportional to i{Iiz,3x — l2z,3x) and, from this, taking Lie brackets with elements in 
C we can obtain all the elements in the basis of Ai indicated in (0). Thus, both C and Ai 
are included in £. A basis of Af can be obtained by Lie brackets of appropriate elements of 
M (possibly adding an element of C). Finally, a basis of TZ can be obtained by Lie brackets 
of appropriate elements of Ai and A/". Therefore the Lie algebra ^ is a subalgebra of C. 
The two Lie algebras coincide if J12 = 0. This is the case remarked above of a system that, 
according to Theorem ^ is state controllable, since C is isomorphic to sp{4), but not operator 
controllable. If Ju 7^ 0, then the matrix A is not in the Lie algebra A. However, it is still 
possible to generate A, which is isomorphic to sp{A) and, applying part 5 of Theorem |I|, we 
conclude that C = su{8) in this case, and the system is operator controllable. 

The results of this section and Section |^ remain true even if we set Uz = constant in the 
model ([l0|) if we assume that there exist no two values for the gyromagnetic ratios 71 and 



72 such that 71 = —72 (cfr. Remark O 



8 Conclusions 

We have presented an analysis of the Lie algebra structure associated to a system of n spin | 
particles with different gyromagnetic ratios and inferred its controllability properties. These 
only depend on the properties of a graph obtained by connecting two nodes representing 
two particles if the coupling constant between the two particles is different from zero. Con- 
trollability of the state and of the unitary evolution operator are equivalent for this class 
of systems. If the system is not controllable then it is a parallel connection of a number of 
controllable systems equal to the connected components of the associated graph. The latter 
result can be easily generalized to the case where the connected components do not represent 
controllable subsystems, which is a case that might occur if some of the gyromagnetic ratios 
are equal. 

We have given a complete description of the low dimensional cases (up to a number of 
particles equal to three) with isotropic interaction and possibly equal gyromagnetic ratios. 
This analysis is of interest since, in many physical situations, a small number of particles is 
controlled. These results also provide an example of a quantum system which is controllable 
in the state but not in its unitary evolution operator. Thus, the equivalence of the two notions 
of controllability, proved for spin systems in the case of different gyromagnetic ratios, is no 
longer true if some of the gyromagnetic ratios are equal. 

This paper also presented a general sufficient condition of controllability for spin systems 
in terms of the associated graph. 
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Appendix A: Some properties of the matrices Ikiii,k2i2,-Air 



We collect in this appendix some properties of the matrices Ikih,k2i2,----,kriri particular 
involving the commutators of these matrices. These relations can be easily proven by using 
the fundamental property: 



where A and B are square matrices of the same dimensions and B and D are square matrices 
of the same dimensions as well. 

Property 1: If {ki, k2, K} fl {^i, ^2, kg} — then 



[A0B,C^D]^ [A, C]^BD + CA^ [B, D], 



(63) 



kih ,k2l2,--;krlr1 ^k\m\,k2' 



>m2,---,ksm, 




(64) 



for every possible combination of Zj's and m^ 's. 

Property 2: Assume that k e {ki, k2, K}, and, in particular, k = kj. 



(a) If Ij — m, then 




(65) 



(b) if [(7i., (7^] = a^T, with a — ±i, then: 



klh,k2l2,---,krlr-l ^ k,m\ J-kih,...,kjT,...,krlr- 



(66) 
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Appendix B: Proof of property ( [55l ) 



In order to see (^51) , we write any element X of u(4) as follows (we use the definition 



o"o := ^2x2 and the ordering < x < y < z) 

X = ^ ajkiaj ® ak= X! ^n'^^i ® + 

j,k=0,x,y,z j=0,x,y,z 
j<k j<k 

From this expression, it is immediate to see that X G JF^ if and only the terms in the last 
sum are all zero. Therefore a basis of J-'p is given by the matrices of the form 

i {ai ®ay + a^® ai) , (67) 



with Z, f = 0, X, y, z. In view of this fact, it is sufficient to verify (|55D on all the matrices of 
the form (|67D. This last fact is only a straightforward calculation. 



Appendix C: Structure of the Lie algebra C in the case 

n = 3, Ji3 = J23 7^ (71 = 72 7^ 73) 

We look at the following vector subspaces of H {H is the vector space defined in (|5^). 



C := span i{Iiv + hvJzw v,w = x,y,z} (68) 

M := span z{/i^,3^ + l2v,3w v,w = x, y, z} (69) 

Af := span i{Iiu,2w,3p + hw,2v,3p v,p,w = x, y, z] (70) 

Q := span i{/i^,2«, + hw,2v v w = x,y,z} (71) 

TZ := span i{Iix,2x - hy,2y, hx,2x - hz,2z} (72) 

The following commutation relations are easily verified 

[C, C] C C, [C, M]CM, [C, AO C A/-, [C, Q] C AT © 7^, [C, TZ] C Q; (73) 

[M,M]CC®M®Af, [M,Af]CM®Q®n, [M,Q]cM [7W,7^]CA^; (74) 

[A/", A/"] C C © AT, [Af,Q]CC®M, [AT, 7^] C (75) 

[Q,Q]CC, [Q,7^]CC; (76) 

[7^,7^] = o. (77) 



A basis of C is generated by the matrices B^., By, B^ according to Lemma [4.1| , while a 
basis in Ai can be obtained calculating the Lie bracket of A with I^^ G C then taking the 
Lie bracket with I^z so as to obtain Iiz,3x + l2z,3x- From this, taking the Lie bracket with 
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elements in C, we can obtain all the elements in the basis of Ai in (|74D . A basis of A/" is 
obtained by Lie brackets of elements in A^. In particular, to obtain elements of the form 
^-?'iD,2«,3t«, we calculate [ili^^3i+il2v,3i, -^hv,3p-ihv,3p]-lihw, "f^ith p ^ lv,p,l e {x,y,z}and 
iaw = [ci, <Jp]. To obtain elements of the form iliv,2w,3p + ihw,2v,3p, v w, v,w,p E {x, y, z}, 
we can calculate Lie brackets of elements of the form ilim,3m + ^-^2m,3m, with elements of 
the form ilin,3n + il2n,3n) with n ^ m and, possibly, calculate the Lie bracket with an 
element of the form 2/3;, l,m,n G {x,y,z}. A basis of Q can be obtained by Lie brackets 
between elements of the form iliv,2v,3x ^ A/" with elements of the form iliw,3x + ihw,3x ^ A/1, 
V ^ w, v,w E {x,y,z}. A basis of 71 can be obtained by the Lie bracket of elements 
ihv,2w,3x + ihw,2v,3x G A/" with elements ilip,3x + il2p,3x, with p v w, p,v,w e {a:, y, z}. 

Notice that if J12 = 0, then A is an element of the Lie algebra above described C (B Ai(B 
Af ® Q®Tl, while if 7^ we have C ® M ® Af ® Q ® TZ = n/span{A}, and the Lie 
algebra C, in this case, coincides with H. 
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